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Abstract 

The soil bacterium Mycobacterium smegmatis is able to scavenge the trace concentrations of H2 present in the atmosphere, 
but the physiological function and importance of this activity is not understood. We have shown that atmospheric H2 
oxidation in this organism depends on two phylogenetically and kinetically distinct high-affinity hydrogenases, Hydl 
(MSI\/lEG_2262-2263) and Hyd2 (I\/1SMEG_2720-2719). In this study, we explored the effect of deleting Hyd2 on cellular 
physiology by comparing the viability, energetics, transcriptomes, and metabolomes of wild-type vs. Aliydl cells. The long- 
term survival of the Aliycl2 mutant was significantly reduced compared to the wild-type. The mutant additionally grew less 
efficiently in a range of conditions, most notably during metabolism of short-chain fatty acids; there was a twofold 
reduction in growth rate and growth yield of the Ahydl strain when acetate served as the sole carbon source. Hydl 
compensated for loss of Hyd2 when cells were grown in a high H2 atmosphere. Analysis of cellular parameters showed that 
Hyd2 was not necessary to generate the membrane potential, maintain intracellular pH homeostasis, or sustain redox 
balance. However, microarray analysis indicated that Ahyd2 cells were starved for reductant and compensated by rewiring 
central metabolism; transcripts encoding proteins responsible for oxidative decarboxylation pathways, the urea cycle, and 
ABC transporter-mediated import were significantly more abundant in the Ahydl mutant. Metabolome profiling 
consistently revealed an increase in intracellular amino acids in the Aliydl mutant. We propose that atmospheric H2 
oxidation has two major roles in mycobacterial cells: to generate reductant during mixotrophic growth and to sustain the 
respiratory chain during dormancy. 
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Introduction 

In recent years, it has emerged that a number of soil 
Actinobacteria of the genera Mycobacterium, Streptomyces, and 
Rhodococcus oxidise the trace concentrations of H2 found in tlie 
lower atmosphere [1,2,3,4]. In addition to being biogeochemically 
important [5], scavenging of tropospheric H2 is physiologically 
unusual; all other characterised hydrogen-oxidising organisms are 
only capable of recycling the high concentrations of H2 evolved 
through other biological processes or geothermal activity [6] . The 
purpose and importance of hydrogen scavenging in the physiology 
oi Actinobacteria nevertheless remains to be understood. It is also 
to be determined whether this process influences the composition 
of microorganisms in soil ecosystems. 

Work in our laboratory has resolved the determinants of 
hydrogen scavenging. The soil bacterium Mycobacterium smegma- 
tis catalyses atmospheric H2 oxidation using two high-affinity, 
membrane-associated, oxygen-dependent [NiFe] -hydrogenases 
[3]. Both of these enzymes are expressed during exponential 
growth, though their expression and activity is significantly higher 
during the transition to stationary phase due to carbon-limitation. 



The fast-acting Group 2a [NiFe] -hydrogenase Hydl 
(MSMEG_2262-2263) is responsible for the majority of whole- 
cell H2 oxidation. In contrast, the Group 5 [NiFe] -hydrogenase 
Hyd2 (MSMEG_2720-2719) is a much slower-acting enzyme in 
whole-cells [7,3]. Despite its low activity, Hyd2 has been shown to 
be important for the growth of M. smegmatis [8] . Furthermore, 
orthologs of this enzyme are more widely distributed among 
sequenced Actinobacteria and are apparently responsible for the 
tropopheric H2 uptake of streptomycetes and rhodococci [9,4]. It 
should also be noted that M. smegmatis also encodes a further 
hydrogenase, Hyd3; this enzyme is only expressed during oxygen- 
limitation, where we propose it serves to couple the reoxidation of 
NAD(P)H to the evolution of hydrogen [7,8]. 

In this work, we provide insight into the physiological role of 
hydrogen scavenging by observing the effect of deleting Hyd2 
throughout exponential growth, upon entry into stationary phase, 
and during long-term survival. Using a combinatorial approach, 
we show that hydrogen scavenging is required for the efficient 
metabolism of certain carbon sources and infer that atmospheric 
H2 is a source of reductant for mycobacterial metabolism. 
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Materials and Methods 

Bacterial strains and growth conditions 

All bacterial strains used in this study are listed in Table SI. 
Mycobacterium smegmatis mc^l55 [10] and derived mutants [7,8] 
were maintained on LB agar plates supplemented with 0.05% (w/ 
v) Tyloxapol (Sigma-Aldrich). For broth culture, M. smegmatis was 
grown in Hartmans de Bont (HdB) minimal medium [1 1] 
supplemented with the stated carbon sources, 0.05% Tyloxapol, 
and 10 (XM NiS04. Cultures were incubated at 37°C with 
agitation (200 rpm) in 30 mL medium in 125 mL aerated conical 
flasks. Culture volumes were upscaled to 500 mL in 2.5 L flasks 
for transcriptome analysis and 100 mL in 500 mL flasks for 
metabolome analysis. Cells were inoculated to an initial optical 
density of 0.005. Optical densities to assess growth were measured 
at 600 nm (ODgoo) in a Jenway 6300 spectrometer. Cultures were 
diluted in 0.85% saline to bring the ODgoo below 0.5 when 
measured in cuvettes of 1 cm light path length. To count colony 
forming units (CFU mL~'), each culture was serially diluted in 
phosphate-buffered saline (PBS) (pH 7.0) and spotted on to agar 
plates [12]. A markerless deletion of the Hyd2 large subunit 
(MSMEG_2719) was complemented with a pOLYG vector 
containing the hyd2 operon (MSMEG_2720-2718) in order to 
minimise disruption to hydrogenase maturation and folding [8] . fi- 
galactosidase assays and amperometric hydrogen measurements 
were performed as previously described [7]. 

Challenge experiments 

For acid challenge experiments, the strains were grown on HdB 
media at pH 7.0 to ODgoo = LO. They were subsequently pelleted 
(7,000 xg, 10 min, RT), washed in 100 mM citrate/phosphate 
buffer (pH 7.0), and resuspended in 100 mM citrate/phosphate 
buffer (pH 3.0 or pH 5.0). All buffer preparations contained 
22 mM glycerol, 0.05% Tween80, and trace metals. Following 
acid challenge, the survival of ceUs was measured by measuring 
colony forming units (CFU mL~'). The minimum inhibitory 
concentrations (MICs) of pH 5.0-challenged cells to the proto- 
nophore carbonyl cyanide m-chlorophenylhydrazone (CCCP) was 
determined using serial dilutions as previously described [13]. 

Measurement of internal pH and membrane potential 

Internal pH and membrane potential was measured in M. 
smegmatis grown on HdB minimal medium at 2 h following the 
induction of stationary phase. Internal pH was calculated by 
determinmg the partitioning of a radioactive probe between 
intracellular and extracellular fractions. Cultures of 1 mL were 
incubated with 11 (iM ['*C] benzoate (10-25 mCi mmol') 
(pH 7.5) (37°C, 10 min) and centrifuged through silicone oil 
(BDH Laboratory Supplies) (16,000 xg, 5 min, RT). A 20 |il 
sample of the supernatant was removed. The tubes were otherwise 
frozen (— 80°C, 60 min) and the cell pellets were removed with 
dog nail clippers. Samples of the supernatant (extracellular 
fraction) and pellet (intracellular fraction) were dissolved in 
scintillation fluid (Amersham). The relative concentrations of 
[ C] benzoate in each sample was measured using a LKB WaUac 
1214 Rackbeta liquid scintillation counter (Perkin Elmer Life 
Sciences). The internal pH was calculated from the uptake of [ C] 
benzoate using the Henderson-Hasselbalch equation as previously 
described [14] . Membrane potential was measured by a equivalent 
method by detemining the partioning of 5 |J,M [' H] methyltri- 
phenylphosponium iodide ([H]TPP^) (30-60 Ci mmol ). The 
membrane potential was calculated from the uptake of [■'H]TPP"'" 
using the Nernst equation [15]. 
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Figure 1. The growth and survival of Mycobacterium smegmatis 
in the presence and absences of hydrogenases. (A) Growth of M. 
smegmatis mc'^155 and hyd mutants into carbon-limitation. (B) Survival 
of M smegmatis mc'^ISS and fiyd mutants during carbon-limitation. The 
strains were grown in aerated conical flasks on HdB minimal medium 
supplemented with 22 mlVI glycerol. Growth is shown in ODsoo- Survival 
is shown in percentage colony forming units relative to day four. 
Legend: Blue circles/bars = Wild-type; Red squares/bars = A/iyc/1 23; 
Orange point-up triangles/bars = Ahyd1; and Purple point-down trian- 
gles/bars = A^yd2. Error bars show standard deviations from biological 
triplicates. * = p<0.05, ** = p<0.01, *** = p<0.001 difference relative to 
wild-type bars (Student's T-test, unpaired, two-tailed). 
doi:1 0.1 371 /journal.pone.01 03034.g001 

IVleasurement of [NAD^]/[NADH] ratios 

1 mL cultures were centrifuged (150,000 xg, 3 min, RT) and 
resuspended in either 0.3 mL 0.2 M HCl (for NAD"^ extraction) or 
0.3 mL 0.2 M NaOH (for NADH extraction). The cultures were 
heated (55°C, 10 min), cooled (0°C, 5 min), and neutralized with 
either 0.3 mL 0.1 M NaOH (for NAD'^ extraction) or 0.3 mL 
0.1 M HCl (for NADH extraction). After centrifugation 
(150,000xg, 3 min RT), the supernatants were coUected. 200 \il 
supernatant was transferred into cuvettes containing 50 |J,1 1 M 
bicine (pH 8.0), 50 10,1 40 mM EDTA (pH 8.0), 50 nl 4.2 mM 3- 
[4,5-dimethylthiazol-2-yl] -2,5-diphenyltetrazolium bromide, and 
50 10,1 16 mM phenazine ethosulfate. NAD* and NADH concen- 
trations were measured by addition of 50 |J,1 ethanol and 5 U yeast 
alcohol dehydrogenase II [16]. The rate of reduction of 3-[4,5- 
dimethylthiazol-2-yl]-2,5-diphenyltetrazoIium bromide was mea- 
sured photometrically at 570 nm and was proportional to the 
concentration of standards of each cofactor. 

RNA extraction 

M. smegmatis mc^l55 and Ahyd2 cells were grown synchro- 
nously in aerated conical flasks. At 1 hour following entry into 
stationary phase due to carbon-limitation, the cells were harvested 
for RNA extraction and microarray analysis. 500 mL of each 
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Figure 2. Observation, complementation, and recovery of mutant growth phenotypes. Strains were grown on HdB minimal medium. (A) 
Growth on 5.5 mM glycerol. (B) Complementation on 5.5 mM glycerol in the presence of 50 |ig mL"' hygromycin. (C) Growth on 12.5 mM acetate. 
(D) Partial complementation on 12.5 mM acetate in the presence of 50 |j,g mL"^ hygromycin. (E) Growth on 12.5 mM acetate in serum vials injected 
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with 10% pure N2. (F) Growth on 12.5 mM acetate in serum vials injected with 10% pure H2. Legend: Blue circles = Wild-type (or wild-type with empty 
pOLYG vector for complementation); Red squares = Ahyd^ 23; Orange point-up triangles = Ahyd^ ; Purple point-down triangles = Ahydl (or Ahydl with 
empty pOLYG vector for complementation); and Grey diamonds = A/iyd2 with pOLYG vector expressing IV1SMEG_2720-2719. Error bars show standard 
deviations from biological triplicates. 
doi:1 0.1 371/journal.pone.01 03034.g002 



culture were mixed with 1000 mL cold glycerol saline (3:2 v/v) (- 
20°C), centrifuged (27000xg, 20 min, -20°C), and resuspended 
in glycerol saline (1:1 v/v) (— 20°C). Cell lysis was achieved by 
three cycles of bead-beating in a Mini-Beadbeater (Biospec) at 
5,000 rpm for 30 sec. Total RNA was extracted using TRIzol 
reagent (Invitrogen) according to the manufacturer's instructions. 
DNA was removed from the RNA preparation by treatment with 
2 U RNase-free DNase using the TURBO DNA-/r«« kit (Ambion), 
according to the manufacturer's instructions. The concentration 
and purity of the RNA was determined using a NanoDrop ND- 
1000 spectrophotometer, and its integrity was confirmed on a 
1.2% agarose gel. 

Microarray analysis 

Transcriptome analysis employed glass slide DNA microarrays 
provided by the Pathogen Functional Genomics Research Center 
(PFGRC), which is funded by the National Institute of Allergy and 
Infectious Diseases. The arrays represented every open reading 
frame of the genome of M. smegmatis mc^l55 with 7,736 unique 
70-mers spotted in triplicate. Samples for microarray analysis were 
prepared and hybridised based on standard operating protocols 
(SOP) M007 and MOOS from The Institute of Genomic Research 
(TIGR) [17]. 5 |a,g extracted total RNA was reverse-transcribed 
and aminoaUyl (aa)-labelled using 3 |J,g random primers (Invitro- 
gen), Superscript III reverse trancriptase (Invitrogen), and a 
25 mM aa-dUTP labelhng mix (2:3 aa-dUTP to dTTP) (Sigma- 
Aldrich). The synthesised cDNA was labelled with cyanine-3 (Cy3) 
or cyanine-5 (Cy5) fluorescent dyes (GE Healthcare Biosciences) 
for 2 h. After measurement of the concentration of cDNA and 
incorporated dyes (NanoDrop ND-1000 spectrophotometer), the 
labelled probes were mixed in equal ratios according to 
instructions in SOP M007. Prior to microarray hybridisation, 
the microarray slides were blocked, washed, and dried as described 
in SOP MOOS. The slides were immediately hybridised with the 
prepared samples and incubated overnight. After hybridisation, 
slides were washed with progressively more stringent buffers and 
dried as per SOP MOOS. Slides were immediately scanned using 
an Axon GenePix4000B microarray scanner (Molecular Devices) 
and analysed with the TM4 suite programs Spotfinder, MIDAS, 
and MeV as previously described [IS]. Gene expression ratio (fold 



change from Ahyd2 vs. wild-type) was calculated from the 
normalised signal intensities. Microarrays were hybridised using 
RNA from each of the four biological replicates. Cy3 and Cy5 dye 
swaps were employed between replicates. 

Quantitative RT-PCR 

cDNA was synthesized from 1 (Xg of RNA for each sample with 
the Superscript III Reverse Transcriptase Kit (Invitrogen). After 
cDNA synthesis, quantitative RT-PCR was performed using 
Platinum SYBR Green qPCR SuperMix-UDG with ROX 
(Invitrogen) according to the manufacturer's instructions. Primers 
(Integrated DNA Technologies) for 10 genes (Table S2) were 
designed with the publicly available Primer3 software. Primer 
pairs were optimised to ensure efficient amplification. The real- 
time PCR reactions were conducted in ABI Prism 7500 (Applied 
Biosystems). Relative gene expression was determined from 
calculated threshold cycle (Cj-) values that were normahsed to 
the gene sigA (MSMEG_275S) as an internal normalisation 
standard. 

Metabolome analysis 

M. smegmatis mc^l55 and Ahyd2 cells were grown synchro- 
nously in aerated conical flasks. At 1 hour following entry into 
stationary phase due to carbon-limitation, samples were collected. 
To prepare samples for extracellular metabolite analysis, 15 mL of 
culture were centrifuged (27000 xg, 10 min, RT) and the 
supernatant was stored at — 20°C. To prepare samples for 
intracellular metabolite analysis, 15 mL of each culture were 
quenched with 30 mL cold glycerol sahne (3:2 v/v) (-20°C) [19], 
centrifuged (27000 xg, 20 min, — 20°C), and resuspended in 1 mL 
glycerol saline (1:1 v/v) (— 20°C) for storage. The samples were 
recentrifuged prior to metabolite extraction and the pellets were 
submitted to metabolite extraction. Before extraction, 20 |J.l of the 
internal standard (10 mM L-alanine-rf4) was added to each 
intracellular and extracellular sample. The metabolites were 
extracted and derivatised as described in existing protocols [20]. 
The intracellular and extracellular metabolites were analysed 
using a gas chromatograph (GC-7S90) coupled to a mass 
spectrometer (MSD5975) (Agilent Technologies) with a quadru- 
pole mass selective detector (EI) operated at 70 eV. The results 



Table 1. Energetic parameters of wild-type and Ahydl cells two hours following the induction of stationary phase. 





Carbon Source 


5.5 mM Glycerol 




1 2.5 mM Acetate 




Strain 


WT 


\hyd2 


WT 


AhycD. 


Growth Yield (ODeoo) 


0.86±0.01 


0.73 ±0.04 


0.37±0.03 


0.19±0.03 


Growth Rate (h^) 


0.19±0.04 


0.17±0.01 


0.17±0.01 


0.07±0.01 


External pH 


6.3±0.1 


6.5±0.1 


7.5±0.1 


7.5±0.1 


Internal pH 


7.3±0.2 


7.1 ±0.1 


7.1 ±0.1 


7.2±0.2 


Membrane Potential (mV) 


-166±7 


-163±5 


-153±4 


-155±9 


NAD^ Concentration (|xM) 


1.6±0.4 


1.6±0.1 


0.90±0.14 


0.95±0.12 


NADH Concentration (|iM| 


3.6±0.5 


2.6±0.5 


0.55±0.14 


0.53±0.12 



Error margins represent standard deviations from three biological replicates, 
doirl 0.1 371 /journal.pone.Ol 03034.t001 
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Figure 3. Genes with significant clianges in expression in t^hydl 
vs. wild-type cells. Both strains were grown synchronously on HdB 
minimal medium supplemented with 22 mM glycerol. RNA was 
extracted at 2 hours following induction of stationary phase. The genes 
were classified as significantly upregulated if expression ratio >2.0, p 
valuesO.05. The genes were classified as significantly downregulated if 
expression ratio <0.5, p valuesO.05. The number of genes affected are 
listed by functional category. The full list of genes in each category is 
shown in Table S4, Table S5, and Dataset S1. 
doi:1 0.1 371 /journal.pone.01 03034.g003 

were processed using R package Metab [21], and samples were 
analysed using the mass spectral deconvolution and identification 
system [20]. The final concentration of metabolites was deter- 
mined using the GC peak intensity of methyl chloroformate 
derivatives. Compounds considered false-positives were eliminat- 
ed, and the intensity of each metabolite was normalised relative to 
the intensity of the internal standard. Samples of five technical 
replicates were collected and analysed from each of the three 
biological replicates. 

Results 

Hydrogen scavenging enhances survival of M smegmatis 

To determine the physiological importance of utilization, the 
growth and survival of Mycobacterium smegmatis mc^l55 and 
hydrogenase mutants was measured in a range of conditions. All 
strains were grown in aerated conical flasks on HdB minimal salts 
medium supplemented with different carbon sources and the non- 
metabolisable detergent Tyloxapol. During growth on 22 mM 
glycerol (pH 7.0), there was no significant difference in the specific 
growth rates (~0. 25 h ) or final growth yields (~ 1 0 CFU mL ) 
of the strain lacking Hydl (Ahydl), Hyd2 {Ahyd2), and aU three 
hydrogenases {Ahydl23) (Figure lA). However, the long-term 
survival of these strains was compromised following the onset of 
carbon-limitation. Each strain lost viability at a significantly and 
reproducibly faster rate than the wild-type. At eight days and at all 
sampling points thereafter, the Aliydl23 strain produced at least 
40% fewer colony forming units compared to the wild-type. For 
example, after 12 days, 3.6 x 10 CFU mL~ were counted for the 
wild-type strain compared to 1.7xlO'C FU mL ~ ' for the Ahyd 1 2 3 
strain The long-term survival of the Aliyd 1 and Ahyd2 strains was 
also reduced (Figure IB). Thus, atmospheric Hj is among the 
substrates that M. smegmatis employs to maintain viability when 
deprived of organic carbon sources. 

M. smegmatis grows optimally by co-metabolising 
organic carbon sources and atmospheric H2 

Major growth phenotypes were observed when the hydrogenase 
mutants were grown at lower carbon concentrations or oxidised 
carbon sources. When the concentrations of glycerol was reduced 
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0.1 1 10 

Relative Abundance 



100 



Figure 4. Profiles of intracellular metabolites in \hyd2 vs. wild- 
type cells. Metabolites were detected by gas chromatography-mass 
spectrometry (GC/MS). The values show the relative abundance of each 
metabolite detected in the samples (arbitrary units) on a logarithmic 
scale. Legend: Green = Wild-type, Yellow = A/iyd2. Means were calculat- 
ed from three biological replicates and five technical replicates for each 
strain, p values were determined using a Student's T-test. * = p<0.05, 
** = p<0.01 difference relative to wild-type bars. 
doi:10.1 371/journal.pone.01 03034.g004 

from 22 mM to 5.5 mM, the Ahydl23 strain grew to a final ODgoo 
approximately 25% lower than the wild-type (Figure 2A). Growth 
of the mutant strain was also defective on a range of other carbon 
sources. There was a twofold reduction in the growth yield of die 
A/i3)(il23 strain compared to the wild-type when the short-chain 
fatty acid acetate was available as the sole carbon source. 
Furthermore, the specific growth rate of the Ahyd\2'5 strain 
(0.03 h~') was fivefold lower than the wild- type (0.17 h~') 
(Figure 2C). When the phenotypes were traced to the single 
mutants, it was revealed that the growth of the Aliydl and Ahyd2 
strains were intermediate to those of the wild-type and Aliyd\23 on 
5.5 mM glycerol (Figure 2A) and 12.5 mM acetate (Figure 2C). 
These phenotypes indicate that these enzymes have overlapping 
roles during growth and survival, but are not redundant. It was 
possible to complement the Ahyd2 phenotypes by expressing the 
genes of the operon MSMEG_2720-2718 from the Hsp60 
promoter-driven hygromycin-resistant shuttle vector pOLYG 
[22] (Figure 2B&D); a reduction of growth rate was stiU observed 
in the complementation strain, likely due to inefficient hydroge- 
nase maturation or differences in gene copy number, but the final 
growth yield was similar to the wild-type strain. 

During growth on 12.5 mM acetate as the sole carbon source, it 
was also possible to partially recover the phenotype of the Ahyd2 
strain, but not the Ahydl or AJiydl23 strains, by growing the cells 
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Figure 5. Altered balance of catabolic/anabolic carbon metab- 
olism in \hyd2 cells. Cells lacking Hyd2 compensate for the loss of 
electrons derived from H2 by increasing oxidation of organic carbon 
sources. There is an increased flux though the tricarboxylic acid cycle 
due to upregulation of enzymes involved in oxidative decarboxylation 
(e.g. ketoglutarate-ferredoxin oxidoreductase) (highlighted in yellow) 
and downregulation of those involved in anaplerosis (i.e. isocitrate 
lyase) (highlighted in blue). We model that loss of CO2 through 
oxidative decarboxylation reactions is principally responsible for the 
decreased biomass of Ahycl2 cells. Oxidative pathways are depicted 
with green arrows, whereas reductive pathways are represented with 
blue arrows. The red text shows the expression ratios of the significantly 
upregulated or downregulated genes in Ahydl vs. wild-type cells. 
doi:10.1371/journal.pone.0103034.g005 

in the presence of 10% (Figure 2E&F). The AIiyd2 strain grew 
at a similar rate to the wild-type (|J,„t = 0.14 h~'; 
M-A/ijd2 - 0.13 h |a,^,jrfi = 0.065 h ') in this condition; these 
cells also attained a higher final growth yield (ODggo wt ~ 0.45; 
ODfioo Ahjd2 = 0.35; ODsoo M;»rfi=0.22) on 10% than when 
ambient H2 was absent (ODeoo wt = 0.40; ODgoo Ahyd2 — 0.2l; 

Ahydl ~ 0.22) (Figure 2F). No differences were observed 
when equivalent volumes of exogeneous N2 were introduced 
(Figure 2E). We model that Hydl (V,^^, (^pp)=12nmol g 
dw ' mill ' [7]) can compensate for the loss of Hyd2 by oxidising 
the majority of this exogenous Hj. However, the converse 
compensation does not occur because the activity of the Hyd2 
(Vmax (app) = 2.5 umol g dw ' min ' [7]) strain is too low to 
consume exogenous H2 at rates meaningful for cell growth. 

We tested whether the type of carbon source had an influence 
on hydrogenase expression or activity. Following growth on 
5.5 mM glycerol compared to 12.5 mM acetate, the three 
hydrogenases were expressed at similar levels (Figure S2A-C) 
and whole-cells oxidised hydrogen at equivalent rates (Figure 
S2D). It was therefore clear that differences in hydrogenase 
activity were not accountable for the differences in the growth 
phenotypes. The Ahyd2 strain grew to proportionately higher 
yields when stationary-phase cells were spiked with 12.5 mM 



acetate. During growth on 50 mM acetate, the Ahyd2 strain also 
remained significantly impaired (|a, = 0.08 h~'; ODgoo fmai- 1-2) 
compared to the wild-type (n = 0.15 h~'; ODgoo fmai— 1-8)- 
Repeated subculturing of acetate-grown cells also caused no 
improvement of growth rates or yields. Thus, the mutant cells 
consumed all the acetate provided in media, but less efficiendy 
coupled its oxidation to growth. 

Ahydl cells maintain redox balance, membrane 
potential, and pH gradients 

It has been shown that both oxidative and evolving hydroge- 
nases can be important for acid tolerance in Enterohacteriaceae 
[23,24,25]. We therefore tested the hypothesis that Hyd2 is 
required for intracellular pH homeostasis in M. smegmatis. 
Consistently, the hydrogenase deletions were significantly im- 
paired compared to the wild-type when the external pH was 
lowered below pH 6.5 (Figure SIA&B); at pH 5.5, there was a 
significant difference in the specific growth rates (|J^ = 0.12 h~'; 
|a,A/,y 123 = 0.05 h~') and final growth yields (ODgoo wt = 0.24; 
ODgoo A/i3if/i23 = 0- 1 1 of the strains (Figure SIC). However, 
the percentage survival, intracellular pH, and protonophore 
susceptibility of exponentially-grown cells challenged at pH 5.0 
or pH 3.0 was similar between the strains (Table S3). It is 
therefore clear that the pH-dependent phenotypes do not reflect a 
defect in pH homeostasis and instead may be a secondary 
consequence of reduced electron input and altered metabolic flux 
in cells unable to scavenge H2. 

Based on previous studies on the physiological roles of uptake 
hydrogenases [6], we predicted that the probable physiological 
role of Hyd2 was therefore to generate a proton-motive force 
through the respiratory chain and/ or reduce coenzymes required 
for reductive processes. To distinguish these possibilities, we 
measured the membrane potential, pH gradient, and NAD^/ 
NADH ratios of wild-type and Ahyd2 strains following growth on 
5.5 mM glycerol or 12.5 mM acetate. However, all three 
parameters were again similar between the wild-type and mutant 
strains and were within normal ranges (Table 1). We hypothe- 
sised that the energetic parameters of the Ahyd2 strain were 
maintained at normal levels due to compensation of any 
deficiencies by organic electron acceptors. 

Organic electron donors compensate for loss of H2 
oxidation 

We performed a microarray analysis to confirm whether 
organic electron donors compensate for loss of tropospheric H2 
oxidation. The transcriptome of Ahyd2 versus wild-type cells was 
compared following growth on HdB minimal medium supple- 
mented with 22 mM glycerol and 0.05% Tyloxapol. These 
conditions were selected because, while Hyd2 is expressed and 
active [7,3], its deletion did not induce a significant growth 
phenotype aside from an extended lag phase (Figure lA); cells 
therefore fully compensate for the loss of the expressed hydrog- 
enases. In this condition, we determined that 65 genes were 
significantiy upregulated (ratio >2.0;^£0.05) (Table S4) and 53 
genes were significantly downregulated (ratio <0.5; ^S0.05) 
(Table S5) in the Aliyd2 strain compared to the wild-type 
(Figure 3; Dataset SI). The genes affected were distinct to those 
previously observed to be affected by slow growth or oxygen- 
limitation [8] . We performed quantitative RT-PCR to confirm the 
quality of the microarray data; the expression ratios of the selected 
genes correlated well with the microarray results (Figure S3). 

The majority of the transcriptional changes involved genes 
implicated in substrate transport and binding, organic acid and 
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alcohol metabolism, and amino acid metabolism (Figure 3). 
There were particularly extensive changes in central intermediary 
metabolism. MSMEG_3706, a bifunctional enzyme that catalyses 
key reactions in the glyoxylate shunt (isocitrate lyase) and 
methylcitrate cycle (methylcitrate lyase) [26], was significandy 
downregulated; this enzyme is usually upregulated during slow 
growth of M. smegmatis and M. bovis BCG [8]. In compensation, 
several predicted glycolytic and tricarboxylic acid cycle enzymes 
were upregulated, i.e. pyrm'ate dehydrogenase, isocitrate dehy- 
drogenase, ketoglutarate-ferredoxin oxidoreductase [27], and 
lactate 2-m()n(K)xygenase [28]. All of these upregulated enzymes 
catalyse oxidative decarboxylation reactions that yield reduced 
cofactors concomitant with the loss of CO2. To compensate for 
downregulation of the methylisocitrate cycle [29], the strain also 
increased expression of the enzymes of the methylmalonyl-CoA 
pathway (propionyl-CoA carboxylase, methylmalonyl-CoA mu- 
tase) that converts propionate to succinate in an ATP-dependent 
manner [30] . These changes suggest that M. smegmatis compen- 
sates for loss of hydrogen oxidation by re-routing carbon flux from 
anabolic to catabolic pathways. 

In amino acid metabolism, the operon encoding the determi- 
nants of the urea cycle was upregulated (MSMEG_3769-3773). 
These ATP-consuming enzymes convert th(- carbon components 
of amino acids into the tricarboxylic acid intermediate fumarate, 
while removing excess nitrogen as urea. Transcripts encoding the 
predicted NAD-dependent glutamate synthase MSMEG_6458- 
6459 were also significantly more abundant. We also observed that 
the putative operons encoding six ABC transporters were 
upregulated, including those predicted to transport trehalose, 
methionine, branched-chain amino acids, and alkane sulfonates. 
Some of these compounds may be scavenged from the cell 
envelope; it has previously been observed that trehalose is 
produced by mycobacteria as a byproduct of mycolic acid cell 
envelope biosynthesis, and the recycling of this compound by a 
homologous ABC transporter is essential for virulence in 
Mycobacterium tuberculosis [31]. 

We observed no significant changes in the expression of any 
hydrogenase-related genes in this condition, including Hydl 
(MSMEG_2263 ratio = 1.5) or Hyd3 (MSMEG_3928 ratio = 0.9), 
suggesting that other hydrogenases do not compensate for loss of 
Hyd2 in this condition. 

Ahydl cells have an altered intracellular metabolome 

The intracellular and extracellular metabolomes of Ahyd2 and 
wild-type cells were determined by gas chromatography-mass 
spectrometry (GC/MS) under the same conditions as the 
microarray. The metabolome profile was generally similar 
between the strains, but there were statistically significant changes 
in the relative abundance of several amino acids and fatty acids 
(Figure 4; Table S6). These changes were consistent with the 
transcriptome data (Table S4; Table S5). The metabolome 
profile substantiates the finding that NAD-dependent glutamate 
synthase is upregulated in Ahyd2 cells. There was a twofold 
increase in the concentration of glutamate in this strain. Consistent 
with the import or recycling of branched-chain amino acids, there 
was also 30-40% more leucine and valine inside Ahyd2 cells 
(Table S6). This is in line with the mutant strains harnessing 
amino acids as electron donors in the absence of Hyd2. 

Discussion 

In conclusion, it is clear that hydrogen scavenging enhances the 
growth and survival oi Mycobacterium smegmatis under a range of 
conditions. Single and double markerless deletions of the 



hydrogen-scavenging enzymes Hydl or Hyd2 grew to lower yields 
than the wild-type strain. Mutant strains were defective when 
ctiltured on minimal medium at low carbon concentrations, acidic 

pH, and, most significantly, on short-chain fatty acids. Reduced 
growth yields of the Ahyd2 strain have also been observed during 
growth on rich media, e.g. LET (lysogeny broth supplemented 
with Tween80) [8]. AU defects were observed when strains were 
grown in flasks aerated with ambient air, i.e. when H2 is available 
at trace concentrations. In high H2 atmospheres, the rapidly- 
oxidising hydrogenase Hydl could compensate for the loss of 
Hyd2, but the Ahydl strain remained dcf(-rti\'(;. M. smegmatis 
therefore grows optimally through mixotrophic metabolism of 
available carbon sources and a ubiqitious hydrogen supply. H2 
scavenging is therefore likely to be a general feature of M. 
smegmatis growth and survival under physiological conditions, and 
may significandy influence the competitiveness of the bacterium in 
soil ecosystems. 

Despite considerable investigation, it has proven challenging to 
resolve why hydrogen metabohsm is so important during 
mixotrophic growth. The membrane potentials, pH gradients, 
and NAD"'"/NADH ratios of the mutant and wild-type strains were 
similar even in phenotype-inducing conditions. This is likely to be 
due to extensive compensation for the loss of hydrogen scavenging 
by organic carbon sources. Transcriptome and metabolome 
analysis inferred that the Ahyd2 strain increased oxidation of 
organic carbon sources through central intermediary pathways. 
There was a pronounced increase in the transcript levels of the 
oxidative decarboxylation reactions of the tricarboxylic acid cycle, 
coupled with a downregulation of the glyoxylate shunt. Increased 
flux through oxidative pathways would increase the production of 
NADH, NADPH, and reduced ferredoxin with concomitant loss 
of CO2 (Figure 5). Mutant cells also upregulated several ABC 
importers for several potential organic electron donors, e.g. 
trehalose, methionine, and branched-chain amino acids, that 
may be scavengc-d from the cell envelope. In addition, metabo- 
lomics ansdysis showed that the amino acids glutamate, valine, and 
leucine were more abundant in the cells of the Ahyd2 mutant. 
Transcriptomics indicate that they may be ultimately broken down 
into fumarate and urea through an upregulated urea cycle. Thus, 
amino acids may also serve as electron donors to compensate for 
loss of hydrogen scavenging. 

To reconcile the available evidence, we propose that tropo- 
spheric H2 principally serves as a source of reductant in carbon 
metabolism during exponential growth of M. smegmatis. The 
electrons yielded from its oxidation may provide the reduced 
compounds required for efficient carbon metabolism. This process 
enables cells to efficiently balance anabolic and catabolic processes 
to maximise yields during mixotrophic growth. Reductant is likely 
to be especially important during growth on acetate. This 
compound is more oxidised than glycerol and its metabolism 
results in higher NAD"'"/NADH ratios (Table 1). In addition, 
acetate molecules cannot be simultaneously used for catabolic and 
anabolic reactions due to loss of CO2. Increased flux through 
catabolic reactions in the Ahyd2 strain on acetate would therefore 
cause reduced biomass via loss of carbon as CO2. Membrane- 
association may be an advantage in this case because it enables the 
hydrogenases to more efEciendy bind the extracellular H2 that 
diffuses into cells. 

We have previously postulated that Hydl and Hyd2 enzymes 
directiy couple the oxidation of H2 to the reduction of O2 via the 
electron transport chain [7]. The membrane association and 
oxygen-dependence of hydrogenase activity indicates it is physi- 
cally and functionally linked to the respiratory chain [3,7]. 
However, it seems improbable that oxidation of nanomolar 
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concentrations of H2 could significantly influence proton-motive 
generation during exponential growth on millimolar concentra- 
tions of carbon sources; our phenotypic and transcriptome studies 
are more consistent with hydrogenases harnessing electrons for 
reductive cellular processes. It nevertheless remains conceivable 
that aerobic hydrogen respiration may be responsible for the 
enhanced long-term survival of wild-type cells compared to 
Ahydl23 cells during carbon-limitation. Hydrogenases are ex- 
pressed at higher levels [8] and oxidise tropospheric H2 more 
rapidly [7] in this condition. Tropospheric H2 oxidation may 
therefore serve as a significant generator of proton-motive force 
when organic carbon supplies are exhausted; Hj is a dependable 
fuel source given it is present at a constant, albeit trace, 
concentration throughout the troposphere. Expression and activity 
profiling suggests that Group [NiFe] -hydrogenases have 
erjui\'alent roles during the sporulation of streptomycetes and the 
adaptation of rhodococci to carbon-limitation [2,4] . 

The processes of using hydrogenases to generate reductant and 
generate proton-motive force need not be mutually exclusive. The 
NADH generated by the Group 3d [NiFe]-hydrogenase of R. 
eulropha, for example, can be simultaneously oxidised in the 
respiratory chain and used as reductant in the Calvin cycle 
[6,32,33]. Tropospheric H2 oxidation may also be coupled to the 
reduction of a multifunctional redox carrier in M. smegmatis. 
Identification and characterisation of the electron acceptors of 
Hydl and Hyd2 is clearly a priority in order to elucidate the 
cellular processes where these enzymes contribute. 

Supporting Information 

Figure SI Importance of hydrogen metabolism for 
growth at acidic pH. AH strains were grown on HdB medium 

supplemented with 22 mM glycerol. pH was adjusted with 
concentrated HCl. (A) Final growth yield of strains grown at a 
range of pHs. (B) Specific growth rate of strains grown at a range 
of pHs. (C) FuU growth curve of strains at pH 5.5. Legend: Blue 
circles = wild-type. Red squares = A^3i(i 123; Orange point-up 
triangles = A/jycJl; Purple point-down triangles = A/i!)i(/2. Error 
bars show standard deviations from biological triplicates. 
(TIF) 

Figure S2 Expression and activity of hydrogenases on 
glycerol and acetate. (A) hydl, (B) hyd2, and (C) hydS 
expression measured with promoter-/flcZ reporter plasmids. M. 
smegmatis mc^l55 harbouring either pJEMhydl-lacZ, p]EMhyd2- 
lacZ, or ]3}EMhyd3-lacZ were grown in HdB minimal medium 
supplemented with either 5.5 mM glycerol or 12.5 mM acetate. 
Samples for [3-galactosidase activity assays were withdrawn in mid- 
exponential (yellow bars) and early stationary (green bars) phase. 
(D) Hydrogenase activity of whole-cells grown in HdB minimal 
medium supplemented with 5.5 mM glycerol or 12.5 mM acetate 
during early stationary phase. Activity was measured amperome- 
tricaUy. Positive values indicate net H2 evolution. Negative values 
indicate net H2 consumption. Error bars show standard deviations 
from three biological replicates. 
(TIF) 

Figure S3 Validation of microarray data by quantitative 
RT-PCR. Microarray data were validated by comparing the gene 
expression changes of selected genes with that of qRT-PCR using 
the same RNA samples. Genes were chosen that were downreg- 
ulated (MSMEG_1203, MSMEG_3706), upregulated 
(MSMEG_3194, MSMEG_3249, MSMEG_3962, 

MSMEG_3769, MSMEG_5059), or unchanged (MSMEG_4640) 
in the microarray. All bars show the expression ratio of genes in 



Ahyd2 vs. wild-type strains. Yellow bars show microarray data. 
Green bars show qRT-PCR data. Error bars represent standard 
deviations from four biological replicates. * —p<0.05, ** =p<O.Ol 
difference relative to zero (Student's T-test). 
(TIF) 

Table SI Bacterial strains and plasmids used in this 
work. 

(DOCX) 

Table S2 qRT-PCR primers used in this study. The 

forward and reverse primers for each gene (MSMEG_XXXX) 

targeted is listed. 

(DOCX) 

Table S3 Intracellular pH homeostasis of M. smegma- 
tis mc^l55 following acid challenge. Percentage survival, 
internal pH, and protonophore susceptibility of wild-type and hyd 

mutants is shown following acid exposure. Cultures were grown on 
HdB supplemented with 22 mM glycerol to OD 1.0. Cells were 
subsequently challenged in 100 mM citrate/pliosphate buffer at 
pH 5.0 or pH 3.0. Error margins show standard deviations from 
three biologically independent replicates. 
(DOCX) 

Table S4 Genes significantly upregulated in \hyd2 vs. 
wild-type microarrays. Means and p values are calculated 
from four microarrays. The genes were classified as significanfly 
upregulated if expression ratio >2.0,p valueSO.05. Less stringent 
criteria was sometimes used when genes were operonic with other 
upregulated genes, or whenp values were perturbed by one clearly 
anomalous replicate. Asterisks are placed next to genes that did 
not meet the strict criteria, but are stiU very likely to be 
upregulated in the Ahyd2 strain. 
(DOCX) 

Table S5 Genes significantly downregulated in ^hyd2 
vs. wild-type microarray. The mean gene expression ratio 
was calculated from tlu' normalised signal intensities of four 
microarrays. The genes were classified as significantiy downreg- 
ulated if expression ratio <0.5 and p value£0.05 (Student's T 
test). Less stringent criteria was sometimes used when genes were 
operonic with other downregulated genes, or when p values were 
perturbed by one clearly anomalous replicate. Asterisks are placed 
next to such genes that did not meet the strict criteria, but are stiU 
very likely to be downregulated in the Ahyd2 strain. 
(DOCX) 

Table S6 List of intracellular and extracellular metab- 
olites in Ahyd2 vs. wild-type cells. Metabolites were detected 

by gas chromatography-mass spectrometry (GC/MS). Values 
show the relative abundance of each metabolite detected in the 
samples (arbitrary units). Means were calculated from three 
biological replicates and five technical replicates for each strain 
and are shown to two significant figures, p values were determined 
using a Student's T-test. Changes were classified as significant 
whenj!)<0.05 (Student's T-test). 
(DOCX) 

Dataset SI All genes with significant changes in expres- 
sion ratios comparing Ahyd2 to WT. See excel file. 
PCLSX) 
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